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/ Description 

Field of the Invention 

s [0001] The present invention relates to determining the position of a mobile station; more specifically, to locating a 
mobile station using time difference of arrival (TDOA). 

Description of the Related Art 

to [0002] A global positioning system (GPS) is commonly used to provide a receiver with accurate measurements of 
its location. The GPS receiver receives a signal from satellites and determines its positions by performing TDOA cal- 
culations based on the known position of the satellites. The receiver is generally attached to a vehicle or boat and is 
provided for this single purpose. The expense of the GPS receivers has generally limited its purchasers to luxury 
vehicle, aircraft, and boat owners. 
is [0003] Digital cellular/PCS phones have become a very convenient and inexpensive way for a person to communicate 
with other persons or communication systems from wherever the person is located. The person can also call 9-1-1 in 
the event of an emergency. However, lo date, wireless communication systems can not accurately determine the lo- 
cation of the caller without the use of satellites and GPS. 

[0004] Current wireless communication systems use multiple access techniques to combine signals from different 
20 sources lo permit many users lo share a common medium without mutual interference. One of the basic lypes of 
multiple access techniques is code division multiple access (CDMA). In CDMA, each base station transmits a pilot 
channel signal, which is essentially an unmodulated pseudo-random noise (PN) sequence. The PN sequence com- 
prises a sequence of PN chips, and each PN chip corresponds to a distance of about 800.4 feet. Each base station 
transmits the pilot channel signal using a different timing offset such that mobile stations can distinguish from which 
2S base station a pilot channel signal was transmitted. 

[0005] The mobile station is time synchronized with a serving base station, i.e., the base station in which the mobile 
station is in communication. The mobile searches time intervals referred to as search windows for the pilot channel 
signals. Each base station is configured to transmit its pilot channel signal such that mobile stations can expect to 
begin receiving no more than one pilot channel signal within each search window. When the mobile station detects a 
30 pilot channel signal, it measures the pilot channel signal strength and records the phase of the pilot channel signal, in 
terms of PN chips, as the pilot channel signal arrives at the mobile station. If the pilot channel signal strength exceeds 
a predetermined threshold, then the base station that transmitted the pilot channel signal is "visible" to the mobile 
station. The measurements and recordings are transmitted from the mobile station to the serving base station or some 
other predetermined location over a reverse link. 
35 [0006] Conventional methods of determining a mobile station's geolocation generally require an indication of dis- 
tances between at least three "visible" base stations and the mobile station. The distance between a base station and 
a mobile station is equal to the time At f for a signal to travel from the base station to the mobile station, multiplied by 
a wave speed u of the signal. If A^u is a distance from the mobile station (having geographic coordinates (Xq.Vq)) to 
a first base station (having known geographic coordinates (x-, ( y-j )), At 2 i> is a distance from the mobile station to a second 
JO base station (having known geographic coordinates (x 2 ,y 2 )), and At distance from the mobile station to a third base 
sta (having known geographic coordinates (x 5 ,y 3 )) I then k on the Pythagorean theorem, the following equations derived 
for a time of arrival (TOA) approach: 
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(1) 



Af 2 V: 



J(x 2 -x Q ) 2 +(y 2 -y 0 ) 2 , 



(2) 



so 



A* 3 v= J(x 3 -x 0 ) 2 +(y 3 -y 0 ) 2 , 



(3) 



to determine the mobile position (Xo,y 0 ). However, in CDMA, the time Atj is unknown because mobile stations have no 
absolute time reference to measure Atj. 
55 [0007] A TDOA approach reduces the number of equations from three to two (equation (3) minus equation (1 ) and 
equation (2) minus equation (1)). The TDOA approach provides accurate location determinations if no system meas- 
urement errors or multi-path effects, described below, are present. Unfortunately, system measurement errors and 
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multi-path effects generally exist and cause deviations from true location determinations. Therefore the above equa- 
tions cannot be used directly to accurately determine the mobile station M's geolocation. 

Summary Of The Invention 

5 

[0008] The present invention addresses these problems by providing a local positioning system (LPS) designed to 
use radio propagation parameters in code-division multiple access (CDMA) forward links or time-division multiple ac- 
cess (TDMA) reverse links to estimate a mobile station's position. 

[0009] The LPS determines the position of the mobile using triangulation methods by minimizing two set of equations, 
10 called cost functions. The first set of cost functions represent distance errors from the "visible" base stations to the 
mobile station, and the second set of cost functions represents position errors in the location estimation of the mobile 
station. Both sets of cost functions include variables common to more than one of the cost functions within the set. 
The cost functions are minimized by estimating values lor the unknown variables within each equation so that the 
distance or position errors in the set are as close to zero as possible, 
rs [001 0] To determine the geographical coordinates of a mobile station when the distances between the mobile station 
and the base stations are not known, the LPS first estimates the distance from the mobile station to the base stations 
to mitigate the system measurement errors and multi-path effect. After the distances are estimated, the LPS estimates 
the geographic coordinates of the mobile station (Xo,y 0 ), based on the estimated distances. 

[001 1] In a preferred embodiment, the LPS is a software implementation on a computer to determine the geographic 
20 location (geolocation) of a mobile station. The LPS receives a data sample includinginformalion indicating arrival limes 
of pilot channel signals at a mobile station and accesses base station information indicating the location of at least 
three cellular or PCS base stations to which the arrival time information is associated. The LPS then estimates the 
distances from the mobile station to the base stations by minimizing a first set of equations or cost functions and 
estimates the geolocation of the mobile station by minimizing a second set of equations or cost functions based on the 
25 estimated distances. 

[0012] The LPS of the present invention provides the benefit of using existing equipment to provide GPS-like posi- 
tioning capabilities. The LPS requires no additional signal detection capabilities, and only requires a minor modification 
to the existing wireless telephone systems. No additional hardware is needed other than the standard CDMA/TDMA 
system, making the LPS cost effective. The LPS can also solve the 9-1-1 mobile location problem for wireless CDMA/ 
30 TDMA systems. Therefore, the LPS can determine the position of a person in distress from their digital phone. 

Brief Description Of The Drawings 

[001 3] The invention will be described in detail with reference to the following drawings, wherein like numerals rep- 
35 resent like elements and: 

Fig. 1 illustrates a mobile station located inside of a triangle formed by three distinct base stations; 

Fig. 2 illustrates a mobile station located outside of a triangle formed by three distinct base stations; 

40 

Fig. 3a is a schematic perspective view of the LPS implementation according to a preferred embodiment of the 
invention; 

Fig. 3b is a schematic perspective view of the LPS implementation according to another preferred embodiment of 
45 the invention; and 

Fig. 4 illustrates a flowchart of a preferred embodiment of the LPS; 

Fig. 5 is a chad illustrating an example performance analysis of the LPS. 

so 

Detailed Description Of The Invention 

[0014] The embodiments described herein arc used in a CDMA forward link triangulation (FLT) system. It is under- 
stood that the embodiments are also applicable to a TDMA reverse link triangulation (RLT) system upon synchronization 
55 of the base stations. 

[001 5] The LPS determines the geographic coordinates of the mobile station by receiving a data sample representing 
information regarding the mobile station, accessing base station information regarding at least three base stations, 
and estimating the location of the mobile station. The LPS determines the location of the mobile station by minimizing 
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a first set of equations or cost functions to estimate the distances between the mobile station and base stations based 
on the data sample and the base station information, and then minimizing a second set of equations or cost functions 
to estimate the geographic coordinates of the mobile station. 

[0016] The LPS is based on TDOA which uses measured phase shift or chip offset information of the pilot channel 

s signals transmitted from particular base stations that are "visible" to the mobile station. A TDOA triangulation approach 
requires time or propagation delay measurements from at least three "visible" base stations. If less than three base 
stations are "visible" to the mobile station, then the LPS will wait for a mobile station report of three "visible" base 
stations or adjust signal strength threshold levels to allow the mobile to recognize more pilot channel signals from other 
base stations. The mobile station frequently measures the pilot channel signal phases so that the location estimation 

10 can be accrued and made more precise over time. 

[0017] Fig. 1 shows a point representing mobile station M located inside a triangle of points representing "visible" 
base stations b 1( b 2 and b 3 at respective distances d 1p d 2 and d 3 from the mobile station M. The distances between 
the base stations are measured as: length b^b 2 between base stations bj and b 2 ; length b 1 b 3 between base stations 
b-, and b 3 ; and length b 2 b 3 between base stations b 2 and b 3 . Angles ot 12 , a 13 and o 23 are formed by arcs ^Mbo, ^Mb^ 

*5 and b 2 Mb 3 respectively. In Fig. 1, angle ct 23 is equal to 360 degrees minus angles ct 12 and a 13 . Fig. 2 is similar to Fig. 
1 except mobile station M is located outside triangle b-,b 2 b 3 and angle a2 3 is equal to angles a 12 plus a 13 . 
[0018] Fig. 3a illustrates a diagram of an LPS implementation. The LPS includes a computer 10 and an article of 
manufacture 20 and may be located at one of the base stations. The article of manufacture 20 includes a computer- 
readable medium and an executable program for locating the mobile station M. 

20 [0019] Fig. 3b illustrates an alternative LPS implementation. The LPS 1 includes the computer 10 for receiving a 
signal 30 carrying the executable program for locating the mobile station M. The signal 30 is transmitted in a digital 
format either with or without a carrier wave. 

[0020] Fig. 4 illustrates a flowchart of LPS for locating the mobile station M in a preferred embodiment. At step S10, 
the LPS 1 reads in data samples (for example, sector number, pilot phase and strength of the pilot channel signal) 

2S from the mobile station M. At step S20, the LPS 1 reads in a cell site tabic which includes information such as the base 
station ID, sector numbers of the base stations, and the base stations' geographic location measured in, for example, 
latitude and longitude. At step S30, the sector numbers of the data samples are matched with those in the cell site 
table to determine from where the pilot channel signals originated. If the pilot channel signals are from at least three 
base stations, then the triangle b 1 b 2 b 3 is formed, as shown in Figs. 1 or 2, and the distances between the mobile station 

30 M and the base stations b y , b 2 and b 3 and the geoiocation of the mobile station M can be determined. 

[0021] The distances between mobile station M and the visible base stations b v b 2 and b 3 are estimated at step 
S40. The computer 10 calculates for distance 6 A such that a set of cost functions for distance errors are minimized 
and determines distances d 2 and d 3 based on the estimated distance d v The estimation of distance d-, and the deter- 
mination of distances d 2 and d 3 based on distance d 1 will be described below. 

35 [0022] The LPS determines the geographic coordinates of mobile station M at S50 using TDOA. The LPS 1 calculates 
the local coordinates of the mobile station M, i.e., (Xq, y 0 ) in relation to the serving base station b 1( and converts the 
local coordinates (X& y 0 ) to global latitude and longitude based on the known latitudes and longitudes of the base 
stations b A , b 2 and b 3 . When succeeding pilot channel signal phase measurements and recordings exist, the geoioca- 
tion of the mobile station M can be re-estimated and averaged to provide an even more accurate analysis. 

40 

Step S40 - Estimating distances between the mobile station and the base stations 

[0023] The two most critical system measurement errors in a TDOA approach are rounding errors in the pilot channel 
signal phase measurement and synchronization errors among base stations. For the pitot channel signal phase meas- 
45 urement, if one chip corresponds to 800.4 feet, then the rounding error (worst case half a chip) contributes to 400.2 
feet in deviation of location. The rounding error can be represented by random variable when satisfying a uniform 
distribution. 

[0024] Ideally, each base station is lime synchronized with the other base stations. Each base station could also be 
time synchronized using a GPS clock. However, the actual clocks in the base stations tend to drift around a nominal 
so value. The drifting error can be represented as a random variable T 2 , satisfying another uniform distribution. The in- 
fluence of the error sources can be added to equal system measurement error T, which is the sum of random variables 
T-, plus T 2 . Accordtngly : a measured pilot channel signal phase pj is equal to a true pilot channel signal phase plus a 
system measurement error T 

[0025] TDOA works best if the measurements being used are those belonging to line-of-sight (LOS) signals because 
55 a straight line is the shortest line between two points. Unfortunately, it is not always possible for the mobile station M 
to receive LOS signals from the base stations b., b 2 and b 3 . A single signal transmitted from any of the base stations 
b v b 2 and b 3 may reflect off different objects such as buildings, trees and vehicles before it reaches the mobile station 
M, and therefore take a longer path than if the signal were a LOS signal. This multi-path effect causes a delay in the 
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arrival of the signal and detrimentally affects the TDOA estimate. 

[0026] Since there is no guarantee that a mobile station M will acquire lineof-sight (LOS) signals from the visible 
base stations b-,. b 2 and b 3 , the delay in arrival time caused by a multi-path signal must be accounted for when using 
TDOA to determine the distance between the mobile station M and the base stations b 1( b 2 and b 3 . However the 

s amount of delay differs depending on the distance and the objects located between the mobile station M and the base 
stations b v b 2 and b 3 and is therefore very difficult to model. Accordingly, a single multi-path parameter u represents 
the proportional time delay caused by the multi-path effect and is modeled as a non-random parameter instead of a 
random number because the single multi-path parameter u must be estimated for all pilot channel signals. Multi-path 
parameter u is generally less than 1, and would equal its maximum of 1 if the mobile station M only acquired LOS 

10 signals from the visible base stations b 1; b 2 and b 3 . 

[0027] It should be noted that one single multi-path parameter u is assumed which implies a homogenous multi-path 
effect. That is, the delay caused by the multi-path effect is assumed to be the same for each pilot channel signal even 
though the multi-path effect on the pilot channel signal from each of the base stations b A , b 2 and b 3 is different. A multi- 
path parameter u, that represent a uniform extra delay can substantially alleviate the multi-path effect. The multi-path 

*s parameter u may be varied in a certain range defined by a model associated with typical environments such as rural, 
urban, suburban, highway, etc. 

[0028] Mobile station M does not know the exact time (as synchronized with GPS) that the base station bj transmits 
a pilot channel signal nor the exact time that the mobile station M receives the pilot channel signal in order to determine 
the time that it takes for the pilot channel signal to travel from the base station bj to the mobile station M. Therefore 

20 ine distances d-,, d 2 and d 3 between the base stations bj and Ihe mobile station M are unknown. 

[0029] However, the base stations are synchronized with each other, and the mobile station M is synchronized with 
the serving base station b v Thus the mobile station M can record chip offsets of pilot channel signal phases emitted 
from remote base stations b 2 and b 3 in relation to the pilot channel signal of the serving base station b v Therefore, 
the mobile M can determine the additional time - after receipt of the pilot channel signal from the serving base station 

2S b-| - required for the pilot channel signals to travel from the remote base stations b 2 and b 3 to the mobile station M 
because the phase of the remote base stations b 2 and b 3 are measurable in relation to the phase of the serving base 
station b h which is set to zero due to the synchronization of the mobile station M with the base station b, . The mobile 
station M identifies a pilot channel signal phase p 2 as a phase difference between the pilot channel signal phase 
recordings of base stations bj and b 2 , and identifies a pilot channel signal phase p 5 as a phase difference between 

30 the pilot channel signal phase recordings of base stations b 1 and b 3 . Accordingly, distance d 2 equals distance d 1 plus 
800.4 feet times the pilot channel signal phase p 2 , or 

d 2 = d 1 +800.4(p 2 )ft (4) 

35 

Similarly, distance d 3 equals distance d-, plus 800.4 feet times the pilot channel signal phase p 3 , or 

d 3 = d 1 +800.4(p 3 )fl (5) 

40 

However, distance d 1 must be estimated before distances d 2 and d 3 can be determined. 

[0030] Consequently, the LPS 1 estimates distance d-, . To search for an estimate for distance 6, , the following equa- 
tions (6)-(8) are cost functions that are minimized for distance errors F 12 , F 13 , and F 23 : 

45 

F 12 = (^1 ^ 2 ) 2 -M 2 ^ 2 -M 2 tf 2 2 +2 M 2 (c/ 1 )(c/ 2 )cosa l2 , (6) 
F 13 = (Z ?1 6 3 ) 2 -M 2 cr i 2 -M 2 d 3 2 + 2 M 2 (o( l )(d 3 )cosa l3l (7) 

so 

F 23 = (d 2 b 3 ) 2 -M 2 d 2 2 -M 2 d 3 2 + 2 M 2 (d 2 )(d 3 )cosa 23 (8) 

for distance d v multi-path parameter u and angles a, and substituting for distances d 2 and d 3 based on equations and 
55 (5). The cost functions for distance errors F 12 , F 1 F 23 must be minimized to arrive at the best 

[0031] The minimization of cost functions F 12 , F 13 , and F 23 can be accomplished using well known minimization 
approaches, such as by steepest decent or incremental search with respect to d v For example, using an incremental 
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search approach, the above cost functions can be minimized by estimating a range for the distance d, and the multi- 
path parameter u, solving the equations (6)-(8) for each predetermined increment in the ranges, and selecting the 
distance d v multi-path parameter u and angles a 12 . cx 13 and a 23 that provide the distance errors F 12 , F 13 and F 23 
closest to zero. After distance d 1 is estimated, distances d 2 and d 3 can be determined using equations (4) and (5). 

5 [0032] Equations (6)-(8) have four unknown values, namely distance 6, , the multi-path parameter u, and angles a 12 
and a 13 . As discussed above, angle is equal to 360 degrees minus angles a 12 and a 13 when the mobile station M 
is located inside triangle b 1 b 2 b 3 . Angle a 23 is equal to angles cx 12 plus a 13 when mobile station M is located outside 
triangle b^b^ However, the angles « 12 and a 13 are determined based on the estimated distance d-,, in other words 
the values of the angles a 12 and oc 13 are determined according the value of distance d v 

w [0033] A skilled practitioner would readily understand that the CDMA (and TDMA) systems can measure a round 
trip delay of a pilot channel signal emitted from the serving base station to the mobile station M and back to the 
serving base station b v This round trip delay provides the benefit of allowing the LPS 1 to use a more narrow scope 
for estimating the range of distance d v 

is Step S50 - Determining the qeolocation of the mobile station 

[0034] After distances d v d 2 and d 3 are estimated, then the mobile station M Cartesian coordinates (Xq ; y 0 ) can be 
estimated by minimizing equations (9)-(11)for cost functions G v G 2 and G 3 : 



20 



25 



G 2 =M 2 (d z y -[(or, -a- 0 ) 2 + (y, -y n yl 



10) 



G 3 = A r(d 3 y -[(.v 3 -,v 0 ) : +(^, - yn y\- 

where Gj, i = 1 , 2, and 3 represents the position error and is zero in an ideal case. However, since distances d v d 2 
35 and d 3 are estimated, equations (7-9) will not be solved exactly, but the best estimate of (x 0 , y 0 ) can be found by 
minimizing Gj. 

Example Estimation and Coordinate Conversion 

40 [0035] The mobile station M is synchronized with the base stations. Consequently, in the mobile station M reply 
message that is sent back to base station b v the phase shift of the reference pilot channel signal transmitted by base 
station b-, is set to zero. The pilot channel signal phases from base stations b 2 and b 3 are recorded in chip off-sets 
from the zero phase shift of base station b v Accordingly, once distance d-, is estimated, distances d 2 and d 3 can be 
determined directly as discussed above. 

45 [0036] In accordance with steps SI 0 and S20 of Fig. 4, the LPS 1 gathers input information including mobile station 
M information and base station b v b 2 and b 3 information. For example, the mobile station M records pilot channel 
signals emitted from base station bj with a base station identifying pilot PN of 432 and a pilot channel signal strength 
of 1 7 (-8.5 dB); from base station b 2 with a base station identifying pilot PN of 76, a pilot channel signal phase p 2 equal 
to 4 PN chips, and a pilot channel signal strength of 2 1 (-1 0.5 dB); and from base station b 3 with a base station identifying 

so pilot PN of 220, a pilot channel signal phase p 2 equal to 3 PN chips, and a pilot channel signal strength of 1 9 (-9.5 dB). 
In accordance with step S30 of Fig. 4, the pilot PNs that are reported by the mobile station M are matched with pilot 
PNs in the sector information stored in a cell site table to determine from which base stations b-j, b 2 and b 3 the pilot 
channel signals were sent. Hero, base station b-, is coll number 138, transmitting a pilot PN of 432 and is located at 
latitude 40.861389 and longitude -73.864167; base station b 2 is cell number 140, transmitting a pilot PN of 76 and is 

55 located at latitude 40.867500 and longitude -73.384722; and base station b 3 is cell number 43, transmitting a pilot PN 
of 220 and is located at latitude 40.878889 and longitude -73.871389. 

[0037] The base station latitudes and longitudes are converted into a local coordinate system (x.y). Base station b-,'s 
coordinates (0.0) are set as the origin, base station b 2 's coordinates (x 2 .0) are set to be on the x-axis, and base station 
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b 3 's coordinates (x 3 , y 3 ) are determined from the known distances among the base stations. 

[0038] In accordance with step S40 of Fig. 4, cost function equations (6)-(8) are then minimized to estimate that 
distance d A =0.801 miles, multi-path parameter n = 0.98. angle a 12 = 1.784084 radians, angle a 13 = 3.002281 radians 
and angle ^ = 1.218859 radians. Based on estimated distance d v distances d 2 and d 3 are determined directly as 

5 described above to equal 0.983620 miles and 0.839603 miles, respectively In accordance with step S50 of Fig. 4, 
equations (9)-(11) are then minimized to determine that the local Cartesian coordinates (Xq y 0 ) equal (0.237018, 
0.357580). These coordinates can be converted back to latitude and longitude so that the mobile station M's location 
can be more easily marked on a map to show which street it is located. In this example, the local Cartesian coordinates 
(0.237018, 0.357580) of the mobile station M's geographic location are converted to latitude 40.867465 and longitude 

io -73.865885. 

[0039] In the above example, angle a 13 is equal to angle a 12 plus angle a^. Therefore the mobile station M is not 
located within triangle b 1 b 2 b 3: but instead is located outside of length b-,b 3 . 

Distance Deviation of Estimates 

75 

[0040] The lower line in Fig. 5 shows an example ot distance deviation (ft) between the true location and LPS esti- 
mated location of mobile station M based on the time deviation (us) caused by system measurement errors, including 
rounding errors in the pilot channel signal phase measurement and synchronization errors. The upper line represents 
the maximum error for performance over a snapshot of time. If a snapshot is extended over time and the distance 
20 deviations are averaged, the distance deviation would become the lower mean error line. Accordingly, if the base 
stations are synchronized, the rounding error in pilot channel signal phase measurement alone approaches 200 feet. 

Reverse Link Trianqulation (RLT) 

25 [0041] In North American TDMA systems, the time arrival is obtained at the base stations rather than at the mobile 
stations. The mobile station transmits a coded digital verification color code (CDVCC) signal as the identity of the 
mobile station. Upon receiving the CDVCC signal the receiving base station time stamps the time of receiving the 
CDVCC signal. If the base stations are synchronized, then the base stations determine the relative time differences 
between arrival of the CDVCC signals by subtracting the time of the receipt of the signal at the first base station from 

30 the time of the later received signals at other base stations. Accordingly, the LPS is applicable to both CDMA and 
TDMA systems. 

[0042] Therefore, equations (6)-(11) can also be applied to TDMA RLT geolocation systems if the clock signals or 
base stations involved in locating a particular mobile station are synchronized. The synchronization could be done by 
installation of GPS. Reverse link signals are transmitted from mobile stations to base stations through the reverse link, . 
35 which is generally a different frequency band than the forward link of CDMA systems, but in a same frequency band- 
and different time slots for TDMA systems. 

[0043] A TDMA reverse link would provide the benefits of better location accuracy if the time arrival is measured at 
the base stations because there would be no chip-rounding error as in a CDMA forward link. In addition, the power 
control in TDMA is not as stringent as in CDMA, therefore making it easier for several base station to "see" signals 
40 from the mobile. Inputs needed by TDMA reverse link triangulation include the identity of the mobile requesting location 
service, the relative time arrivals at the base stations, location (latitude/longitude) of all base stations, and the roundtrip 
delay (measured continuously in TDMA for time alignment purposes). The strength of signal from the mobile station 
is also desired and can be measured at neighboring base stations for assistance of handoff. 

[0044] While this invention has been described in conjunction with specific embodiments thereof, it is evident that 
45 many alternatives, modifications, and variations will be apparent to those skilled in the art. Accordingly, the preferred 
embodiments of the invention as set forth herein are intended to be illustrative, not limiting. Various changes may be 
made without departing from the spirit and scope of the invention as defined in the following claims. 



50 Claims 

1. A method for determining the position of a mobile station, comprising the steps of: 

(a) receiving pilot channel signal information indicating arrival times of pilot channel signals at the mobile 
55 station; and 

(b) estimating the location of the mobile station by minimizing a set of position error cost functions based on 
the pilot channel signal information and on base station information indicating the location of a plurality of base 
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stations to which the arrival times are associated. 
2. The method of claim 1 , wherein the position error cost functions are derived from the equations: 



G x =M 2 {d x y - [(x x -x Q ) 2 +<>, -y 0 ) 2 \, 



G 3 =/J 2 (d 3 ) 2 -[(x 3 -x Q ) 2 +(y 3 -y 0 ) 2 ] 

where p. is a multi-path effect parameter, d 1 is a distance from the mobile station to the first base station, d2 is a 
20 distance from the mobile station to the second base station, d 3 is a distance from the mobile station to the third 

base station, (x 0 ,y 0 ) are the local Cartesian coordinates of the mobile station, (x^y-,) are the local Cartesian coor- 
dinates of the first base station, (x 2 ,y 2 ) are the local Cartesian coordinates of the second mobile station, and (x 3 , 
y 3 ) are the local Cartesian coordinates of the third mobile station. 

25 3. The method of claim 2, wherein 

d 2 = d, + 800.4(p 2 )ft : where p 2 is a phase difference between the pilot channel signal phase recordings of the 
first and second base stations, and 

30 d 3 = d 1 + 800.4(p 3 )ft, where p 3 is a phase difference between the pilot channel signal phase recordings of the 

first and third base stations. 

4. The method of claim 1 , prior to step (b), further comprising the steps of: 

35 (c) receiving the base station information indicating the location of the plurality of base stations; and 

(d) matching the pilot channel signal information to the base station information based on a source identifier 
common to both the pilot channel signal information and the base station information. 

40 5. The method of claim 1 , wherein the arrival times correspond to synchronized timing of the base stations. 

6. The method of claim 1 , prior to step (b), further comprising the steps of: 

(e) estimating a distance from the mobile station to one of the base stations by minimizing a set of distance 
error cost functions including angles formed by the base stations and the mobile station. 



7. The method of claim 6, wherein the distance error cost functions are derived from the equations: 

F 12 =(£? 1 £? 2 ) 2 -|i 2 d 1 2 -u 2 d 2 2 +2u 2 (o' 1 )(or 2 )cosa 12! 



F 23 ^(b 2 b 3 f-\x 2 d 2 -\x Z d 2 +2\x 2 {d 2 ){d 3 )cosa 23 

where b-,b 2 is a distance from the first base station to the second base station, b 1 b 3 is a distance from the first 
base station to the third base station, b 2 b 3 is a distance from the second base station to the third base station, u 
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is a multi-path effect parameter, is a distance from the mobile station to the first base station, d 2 is a distance 
from the mobile station to the second base station. d 3 is a distance from the mobile station to the third base station, 
a 12 denotes the angle formed by lines between the mobile station and a first and second base stations. a 13 denotes 
for the angle formed by lines between the mobile station and the first and a third base stations, and 003 denotes 
5 the angle formed lines between the mobile station and the second and third base stations. 

8. The method of claim 1 , wherein the pilot channel signal information includes at least one of a source identifier, a 
pilot channel signal phase and a pilot strength. 

10 9. The method of claim 8, wherein the base station information includes at least one of the source identifier and base 
station location. 

10. The method of claim 1, further comprising the step of averaging the estimation of the mobile station location with 
a previous estimation of the mobile station to acquire a mean estimation of mobile station location. 

75 

11. The method of claim 1 , wherein the pilot channel signal information is included in a CDMA signal. 

12. The method of claim 1 , wherein the pilot channel signal information is included in a TDMA signal. 

20 13. A local positioning system implemented on a computer for determining the position of a mobile station, comprising: 

means for receiving pilot channel signal information indicating arrival times of pilot channel signals at the 
mobile station; and 

25 means for estimating the location of the mobile station by minimizing a set of position error cost functions 

based on the pilot channel signal information and on base station information indicating the location of a 
plurality of base stations to which the arrival times are associated. 



30 



35 



40 



45 



14. The local positioning system of claim 1 3, wherein the position error cost functions are derived from the equations: 

g, = /r (</,)' -[(*, - vr +(>, -yS'\ 

G 2 =m 2 W,Y -fo -*o) 2 +(y 3 ->'o) 2 L 
C, =M 2 (cf,) 2 -[o 3 -a- 0 ) 2 +0> 3 -y 0 ) 2 ] 



where u is a multi-path effect parameter, d 1 is a distance from the mobile station to the first base station, is a 
distance from the mobile station to the second base station, d 3 is a distance from the mobile station to the third 
base station, (x 0 ,y 0 ) are the local Cartesian coordinates of the mobile station, (x v y 1 ) are the local Cartesian coor- 
dinates of the first base station, (x 2t y 2 ) are the local Cartesian coordinates of the second mobile station, and (x 3 , 
so y 3 ) are the local Cartesian coordinates of the third mobile station. 

15. The local positioning system of claim 14, wherein 

d 2 = d 1 + 800.4(p 2 )ft, where p 2 is a phase difference between the pilot channel signal phase recordings of the 
55 first and second base stations, and 

d 3 - d, + 800.4(p 3 )ft, where p 3 is a phase difference between the pilot channel signal phase recordings of the 
first and third base stations. 
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16. The local positioning system of claim 13 3 prior to the means for estimating the location of the mobile station, further 
comprising: 

means for receiving the base station information indicating the location of the plurality of base stations; and 

■ s 

means for matching the pilot channel signal information to the base station information based on a source 
identifier common to both the pilot channel signal information and the base station information. 

17. The local positioning system of claim 1 3, wherein the arrival times correspond to synchronized timing of the base 
10 stations. 

1 8. The local positioning system of claim 1 3 : prior to the means for estimating the location of the mobile station, further 
comprising: 

means for estimating a distance from the mobile station to one of the base stations by minimizing a set of 
is distance error cost functions including angles formed by the base stations and the mobile station. 

19. The local positioning system of claim 18, wherein the distance error cost functions are derived from the equations: 

F 13 =(t> 1 b 3 f'[i 2 d^ 2 ^ 2 d 3 2 +2^ 2 (d A )(d 3 )cosa A31 

F 23 ^(b 2 b 3 f-\i 2 d 2 2 -^ 2 d 3 2 +2^ 2 (d 2 )(d 3 )cosa 2Z 

where b-,b 2 is a distance from the first base station to the second base station, b-,b 3 is a distance from the first 
base station to the third base station, b 2 b 3 is a distance from the second base station to the third base station, u 
30 is a multi-path effect parameter, d-, is a distance from the mobile station to the first base station, d 2 is a distance 

from the mobile station to the second base station, d 3 is a distance from the mobile station to the third base station, 
a 12 denotes the angle formed by lines between the mobile station and a first and second base stations, a, 3 denotes 
for the angle formed by lines between the mobile station and the first and a third base stations, and 033 denotes 
the angle formed lines between the mobile station and the second and third base stations. 

35 

20. The local positioning system of claim 13, wherein the pilot channel signal information includes at least one of a 
source identifier a pilot channel signal phase and a pilot strength. 

21. The local positioning system of claim 20, wherein the base station information includes at least one of the source 
40 identifier and base station location. 

22. The local positioning system of claim 13, further comprising means for averaging the estimation of the mobile 
station location with a previous estimation of the mobile station to acquire a mean estimation of mobile station 
location. 

45 

23. The local positioning system of claim 1 3, wherein the pilot channel signal information is included in a CDMA signal. 

24. The local positioning system of claim 1 3, wherein the pilot channel signal information is included in a TDMA signal. 

so 25. An executable program embodied on a computer-readable medium for determining the position of a mobile station 
comprising: 

a receiving source code segment for receiving pilot channel signal information indicating arrival times of pilot 
channel signals at the mobile station; and 

55 

an estimating source code segment for estimating the location of the mobile station by minimizing a set of 
position error cost functions based on the pilot channel signal information and on base station information 
indicating the location of a plurality of base stations to which the arrival times are associated. 
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26. The executable program of claim 25, further comprising: 

a computing source code segment for estimating a distance from the mobile station to one of the base stations 
by minimizing a set of distance error cost functions including angles formed by the base stations and the mobile 
station. 



27. A computer data signal comprising: 

a receiving signal segment for receiving pilot channel signal information indicating arrival times of pilot channel 
signals at the mobile station; and 



an estimating signal segment for estimating the location of the mobile station by minimizing a set of position 
error cost functions based on the pilot channel signal information and on base station information indicating 
the location of a plurality of base stations to which the arrival times are associated. 



*5 28. The computer data signal of claim 27 ; further comprising: 

a computing signal segment for estimating a distance from the mobile station to one of the base stations by 
minimizing a set of distance error cost functions including angles formed by the base stations and the mobile station. 

29. The computer data signal of claim 27, wherein the computer data signal is embodied on a carrier wave. 
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